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ABSTRACT 

Mass estimates of K giants are generally very uncertain. Traditionally, stellar masses of single field 
stars are determined by comparing their location in the Hertzsprung-Russell diagram with stellar 
evolutionary models. Applying an additional method to determine the mass is therefore of significant 
interest for understanding stellar evolution. We present the time series analysis of 11 K giants recently 
observed with the WIRE satellite. With this comprehensive sample, we report the first confirmation 
that the characteristic acoustic frequency, i^max, can be predicted for K giants by scaling from the 
solar acoustic cut-off frequency. We are further able to utilize our measurements of t'max to determine 
an asteroseismic mass for each star with a lower uncertainty compared to the traditional method, for 
most stars in our sample. This indicates good prospects for the application of our method on the vast 
amounts of data that will soon come from the COROT and Kepler space missions. 
Subject headings: stars: fundamental parameters — stars: oscillations — stars: interiors 
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1. INTRODUCTION 

According to theoretical calculations of stellar evolu- 
tion, stars of a large mass range, corresponding to main- 
sequence spectral classes from late B to K, all end up 
in roughly the same part of the Hertzsprung-Russell (H- 
R) diagram when they evolve to become red giants with 
spectral classes from late G to M. Here, the evolution 
tracks are closely spaced, which makes it difficult to es- 
timate the mass, and hence to determine the progenitor 
of a red giant star from its position in the H-R diagram. 

The investigation of stellar oscillations (asteroseis- 
mology), in particular solar- like oscillations, provides a 
unique tool to probe the interior and hence also the mass 
of stars. Much excitement therefore followed the first 
clear evidence of solar-like oscillations in a red giant star 
(Frandsen et al. 2002). Subsequently, researchers have 
seen evidence of oscillations in about a dozen K and late 
G type giants, both in nearby field stars (De Ridder et al. 
2006; Barban et al. 2007; Tarrant et al. 2007) and mem- 
bers of the open cluster M67 (Stello et al. 2007). 

However, it remains uncertain if accurate mode fre- 
quencies are possible to obtain, since observational indi- 
cations are that mode lifetimes could be too short (Stello 
et al. 2006), in disagreement with theory (Houdek & 
Cough 2002). Further, these stars might pulsate pre- 
dominantly in radial mode overtones as suggested by 
Christensen-Dalsgaard (2004) or indeed show additional 
non-radial pulsations as indicated by Hekker et al. (2006) 
and Kallinger et al. (2007). Without a full understand- 
ing of which modes we observe, we cannot exploit the 
full potential of the asteroseismic analysis. 

In this paper we aim at utilizing the characteristic 
acoustic frequency, denoted t'max- This asteroseismic 
quantity is relatively insensitive to the mode lifetime and 
the number of excited modes, but still enables us to ob- 
tain information about the stellar interiors. In particular, 
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we will use t'max to estimate the masses of a sample of 
K giants observed with the star tracker on the Wide-Field 
Infrared Explorer ( WIRE) satellite. 

2. OBSERVATIONS AND DATA REDUCTION 

During its lifetime WIRE observed more than 40 
evolved stars, of which we have selected a subset of 11 
K giants that all have long time series and good noise 
characteristics. The time series span from 15 to 61 days, 
and the observations were obtained from February 2004 
to May 2006. In Fig. 1 we show the location of the giants 
in the H-R diagram, together with stars that have previ- 
ously been reported to show evidence of solar-like oscilla- 
tions. We extracted the complete set of evolution tracks 
(M = 0.5-10 M0) from the BaSTI database (Pietrinferni 
et al. 2004), based on their alpha-enhanced standard so- 
lar models without overshooting (Z=0.0198, Y=0.273). 
We only plot a representative subset of tracks in Fig. 1. 

We obtain the raw light curves using the pipeline de- 
scribed by Bruntt et al. (2005). Data reduction proceeds 
in two stages. First, we remove obviously aberrant points 
from the time series. This is done by examination of 
instrumental magnitude, FWHM of the stellar profile, 
centroid position, and background level as a function of 
time. Generally, the majority of data points examined in 
this way fall into a well-defined group with a small per- 
centage of outliers, which are removed. We next remove 
effects from high levels of scattered light at the beginning 
of each orbit by phasing each time series at the satellite 
orbital period, and then subtract a smoothed version of 
the phased light curve. The resulting mean-subtracted 
time series is used for the analysis described below. 

The WIRE star tracker obtained images of each star 
with a cadence of 2 Hz, providing a few million observa- 
tions per star. However, with the expected long oscilla- 
tion period of the K giants (P > 6 hours) we are able to 
speed up the Fourier analysis (see Sect. 4) by binning the 
data in intervals equal to the satellite orbital period (~90 
minutes), corresponding to 230-830 data points per star. 
The resulting time series have a continuous sampling of 
one data point per orbit with only very few small gaps. 
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which translates into a Nyquist frequency of 90 /uHz. 

3. STELLAR PARAMETERS 

To facilitate our investigation we have derived the stel- 
lar parameters for each star. This enables us to predict 
the characteristic acoustic frequency, fmax,pre, where we 
would expect to see excess power in the Fourier spec- 
trum. The stellar parameters are listed in Table 1, and 
sorted according to t'max.obs) which is the frequency of 
the observed excess power (see Sect. 4). 

In the following we explain each column of Ta- 
ble 1. The V magnitude is obtained from the SIMBAD 
database, derived as the mean of the listed values from 
up to six sources, and the standard deviation is adopted 
as a conservative uncertainty. The infrared K band mag- 
nitude is taken from 2MASS (Cutri et al. 2003), the par- 
allax, IT, is from the new Hipparcos release (van Leeuwen 
2007), and the effective temperature is derived using the 
Tcti-{V - K) relation by Alonso et al. (1999). The in- 
ternal error of the T'cfr-(V" — K) relation is only 25 K. 
However, we adopt cTT^g = 100 K as a realistic uncer- 
tainty, in agreement with Kucinskas et al. (2005). The 
color-temperature relation requires as input [Fe/H], but 
is insensitive to log 5. For five stars in our sample, spec- 
troscopic information is available within a uniform col- 
lection (Mc William 1990), and they all agree with hav- 
ing solar metallicity. Hence, we assume solar metallicity 
for all our targets, which is reasonable for such nearby 
stars. The spectroscopically determined effective tem- 
peratures are in good agreement with those derived in 
Table 1. Three stars have interferometrically calibrated 
Tefi (Blackwell & Lynas-Gray 1998; di Benedetto 1998), 
which also agree with our quoted values. 

We then derive the luminosity from log(L/L0) = 
-[V + BC- 51og(1000/7r) + 5 - Mboi,0]/2.5, where BC is 
the bolometric correction obtained from the BC-{V — K) 
relation by Alonso et al. (1999), and Mboi,© = 4.75 is 
the solar absolute bolometric magnitude (recommenda- 
tion of lAU 1999). Since these are all nearby stars we 
neglect interstellar absorption, Ay, but we do include 
ill ths error budget. In general, the relative con- 
tributions to ctl from the individual micertainties are 
o'BC > c^Av > cttt > cy, with asc being dominated by 
<tt^„, while (Ttt and cry are negligible for most stars. 

The quoted range in the photometric mass, A{fphot(see 
Table 1), corresponds to the lowest and highest masses 
of all tracks that go through the la-error box in the H-R 
diagram, while taking mass loss into account in the stel- 
lar models (mass loss parameter rj = 0.4; see Pietrinferni 
ct al. (2004) and references herein). We note that this 
approach will always underestimate the true mass range 
within the Icr-error box due to the non-zero mass step, 
AM, in our grid of tracks, which is AM = 0.1 Mq for 
0.9 < M/Mq < 2.0, AM = 0.2 Mq for 2.0 < M/M© < 
3.0, and AM = 0.5 M© for 3.0 < M/M© < 10.0. Hence, 
if a track is just outside the error box, our quoted mass 
range will be smaller by almost one mass step. Within 
the age of the universe, none of the tracks with progen- 
itor mass M < 0.85 M© and solar metallicity have yet 
evolved to the red giant phase (sec Fig. 1). Hence, we use 
the 0.90 M© track to define the minimum masses'^. It is 

^ We note that for sub-solar metallicities we would expect to see 
red giants with masses below 0.85 Mq. 



beyond the scope of this Letter to go into detail about ad- 
ditional systematic errors in these mass estimates, origi- 
nating from metallicity, overshooting, and differences in 
evolution codes. 

Finally, we calculate the frequency, t'max.pre (see Ta- 
ble 1), where the highest excess power from solar- like 
oscillations is expected in the Fourier spectrum. This is 
obtained by scaling the acoustic cut-off frequency of the 
Sun (Brown et al. 1991) 

_ M/M© (Teff/5777K)3-5 

i^max — J- /J '^max,© , 

where we use the solar value, fc'max,© = 3021 ± 27 /xHz, 
found from 10 independent 30-day time series from the 
VIRGO instrument on board the SOHO spacecraft, us- 
ing the same approach as for the K giants (see Sect. 4). 
Equation 1 has been shown to give very good estimates of 
the frequency of maximum oscillation power based on ob- 
servations of mostly less evolved stars with relatively well 
constrained masses (Bedding & Kjeldsen 2003) compared 
to our targets. As with Mphot, we also quote a range for 
t^inax,pre, which statcs thc extreme values of fmax.pic that 
are within the la-error box in the H-R diagram. To il- 
lustrate this we show an H-R diagram close-up of our 
target stars in Fig. 2. The black dots correspond to the 
stellar parameters given in Table 1 (heading: "Derived" ), 
and the la-error boxes and HD numbers arc also shown. 
We plot L, Teff for selected values of fmax,pre along each 
evolution track, which clearly shows the complexity of 
estimating t'max from the stellar position in this part of 
the H-R diagram. Note that, similar to Fig. 1, only a se- 
lected sample of all evolution tracks are plotted. Along 
each evolutionary track, a given value of t'max,pre (shown 
as identical symbols) can occur up to three times, twice 
on the red giant branch (ascending and descending), and 
once while ascending the asymptotic giant branch. In 
any region where the tracks cross, the mass, and hence 
i^max, is not uniquely determined by the location in the 
H-R diagram. The last two columns in Table 1 will be 
explained in the following section. 

4. ASTEROSEISMIC ANALYSIS 

To look for evidence of solar-like oscillations, we first 
calculate the Fourier spectrum of each star, which is 
shown in amplitude (=^power) in Fig. 3. The mono- 
tonic green curve is a fit to the noise, described by 
a(i^) = a/v'^ + a^jj in power. The parameters a and a^j^ 
are determined following the approach by Stello et al. 
(2007). In addition, we smooth each spectrum (red 
curve) to remove the detailed structure of the excess 
power. Smoothing was done using a moving box av- 
erage twice with a width equal to twice the expected 
frequency separation of adjacent radial modes derived 
as: Au = (M/M©)0-^(i?/R©)-i-5l34.92 /xHz (Kjeldsen 
& Bedding 1995). We note that the location of the ex- 
cess hump in the smoothed spectrum does not depend 
strongly on the mass adopted in the calculation of Ai^. 

From Fig. 3 we see a clear trend of excess power shift- 
ing to higher frequency from top to bottom. A power 
excess at higher frequencies generally corresponds to less 
luminous stars, but this trend is modulated by mass and 
temperature, in agreement with Eq. 1. The gray shaded 
areas in each panel indicate the predicted region of the 
excess power (Table 1, G0I.8). We further note that the 
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amplitude seems to decrease as the width of the envelope 
increases. Similar trends can be seen in Kjeldsen et al. 
(2005) for less luminous stars. 

Despite the large uncertainty in the stellar masses, 
-^photi and hence in t'max.pro our results in Fig. 3 con- 
firm that Eq. 1 is valid for K giants in a large luminosity 
range. Now, if wc turn the argument around, assuming 
that Eq. 1 is exact, and hence interpret any observed 
deviation from this relation to be largely due to an in- 
accurate "photometric" mass, we can use it to infer an 
"asteroseismic" mass. To obtain the asteroseismic mass 
we measure Vrnax by first subtracting the noise fit from 
the smoothed spectra, and then we locate the maximum 
of the residual (see Table 1, Col. 9). By smoothing we 
obtain a more robust measure than trying to locate the 
strongest oscillation mode. To estimate the uncertainty 
of i^maxjobs, wc make 20 simiilations of each time series 
following the approach by Stello et al. (2004). The simu- 
lations show that t'max.obs has an uncertainty of roughly 
10%, but it varies somewhat from star to star, due to dif- 
ferences in pulsation characteristics, the duration of the 
time series, and the noise level. We found that our results 
do not depend strongly on the adopted mode lifetime, r, 
in the range 3d < r < 20d, that wc investigated. We 
then finally determine Mgeis using Eq. 1. 

Our results show that for stars located in the region 
where evolution tracks cross, Alsds has a lower la un- 
certainty than the Mphot mass range (assuming Eq. 1 is 
exact). In regions without crossings the benefit from hav- 
ing measured fmax,obs is less obvious. For those stars, we 
need longer time series to obtain a lower uncertainty in 
t'max,obs, which generally dominates the uncertainty on 
our mass estimate. For a few stars we arc able to mea- 
siu-c the large frequency separation, Au, which poten- 
tially can give a more precise mass estimate than i/max- 
However, for our present data i/max provides the smallest 
mass uncertainties. 

5. CONCLUSIONS 

We have analyzed photometric time series of 11 nearby 
K giants obtained with the WIRE satellite. The Fourier 
transforms show clear evidence that oscillation power 
shifts to higher frequencies for less luminous stars, as an- 
ticipated from scaling the solar frequency of max;imum 
power, z^max,©- We were able to measure fmax and made 
simulations of each star to obtain a realistic uncertainty 
of the measurement. Using a simple scaling relation, 
which relates this frequency to the stellar parameters 
Teff, L/Lq, and M/Mq, we estimated an asteroseismic 
mass. For several stars this approach provides a signifi- 
cantly lower uncertainty of the mass relative to the clas- 
sical mass estimate based purely on comparing stellar 
evolution tracks with the location in the H-R diagram. 

These results show exciting prospects for the cur- 
rent COROT mission (Baghn & The COROT Team 
1998) and the upcoming Kepler satellite (Christensen- 
Dalsgaard et al. 2007), which will both provide much 
more extended times series than WIRE. With longer time 
series we can potentially acqiiire lower uncertainties in 
the fmax measurements, and hence more precise mass es- 
timates. This will indeed be possible with Kepler, which 
will obtain parallaxes, and hence luminosities, of its tar- 
get stars. Our approach could be particularly valuable 
in cases where the oscillation power does not allow accu- 



rate detection of the large frequency separation, Av, due 
to an insufficient number of modes with high signal-to- 
noise. This might, in fact, include most faint solar-like 
pulsators, as well as bright giant stars if indeed the mode 
lifetime does not increase with increasing oscillation pe- 
riods (Stello et al. 2006). 
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thank Tim Bedding to comments on the manuscript. 
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TABLE 1 

Stellar parameters of the WIRE K giants. 
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*The adopted Icr uncertainty of T^g is 100 K. 

''Numbers in parentheses are uncertainties, e.g. for HD41047 the V magnitude and its uncertainty is 5.543 ± 0.005 mag. 
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Fig. 1. — H-R diagram with our 11 target stars including Icr-error boxes. Additional bright field stars (filled star symbols) and M67 
cluster members (empty star symbols) that show evidence of solar-like oscillations axe marked. Dashed lines indicate the approximate 
location of the classical instability strip. Solid lines are evolution tracks. 
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Fig. 2. — H-R diagram of region around our target stars (solid dots inside Icr-error boxes). Gray lines are a representative sample of 
evolution tracks. The additional symbols show points along the evolution tracks corresponding to a given value of i^max.pre (values are 
indicated for each symbol). 
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Fig. 3. — Amplitude spectra of the 11 K giants (HD numbers are indicated). Note the increasing amplitude scale on the ordinate. In 
each panel the red solid line is a smoothed version of the spectrum and the green monotonically decreasing line is a fit to the noise. Gray 
shaded areas indicate the frequency interval of Vmax.prc. The dotted line marks one third of the WIRE orbital frequency, and CTwn is the 
square root of the mean power in the range: 70-80 /xHz. The inset is the spectral window in amplitude, shown in the same horizontal scale 
as the main panels. 



